DBI Galileons in the Einstein Frame: Local Gravity and Cosmology 



Miguel Zumalacarregui^''^, Tomi S. Koivisto^, and David F. Mota^ 
^ Instituto de Fisica Teorica IFT-UAM-CSIC, Universidad Autonoma de Madrid, 
C/ Nicolas Cabrera 13-15, Cantoblanco, 28049 Madrid, Spain and 
^ Institute for Theoretical Astrophysics, University of Oslo, P.O. Box 1029 Bhndern, N-0315 Oslo, Norway 

It is shown that a disformally coupled theory in which the gravitational sector has the Einstein- 
Hilbert form is equivalent to a quartic DBI Galileon Lagrangian, possessing non-linear higher deriva- 
tive interactions, and hence allowing for the Vainshtein effect. This Einstein Frame description 
considerably simplifies the dynamical equations and highlights the role of the different terms. The 
study of highly dense, non-relativistic environments within this description unravels the existence of 
a disformal screening mechanism, which represents an alternative way to investigate the Vainshtein 
effect. Disformal couplings to matter also allow the construction of Dark Energy models, which 
behave differently than conformally coupled ones and introduce new effects on the growth of Large 
Scale Structure over cosmological scales, on which the scalar force is not screened. We consider a 
simple Disformally Coupled Dark Matter model in detail, in which standard model particles fol- 
low geodesies of the gravitational metric and only Dark Matter is affected by the disformal scalar 
field. This particular model is not compatible with observations in the linearly perturbed regime. 
Nonetheless, disformally coupled theories offer enough freedom to construct realistic cosmological 
scenarios, which can be distinguished from the standard model through characteristic signatures. 
The use of more general disformal transformations provides even further relations between scalar- 
tensor theories of gravity. 

PACS numbers: 95.36.-l-x, 04.50. Kd, 98.80.-k 



I. INTRODUCTION 

In the standard ACDM model of cosmology J], the 
universe at the present day appears to be extremely fine 
tuned. The energy scale of the A component is extremely 
small compared to the nai've quantum corrections [5J [3] , 
and yet large enough to be detectable through its effect 
on the cosmological expansion [4H6|. This mystery has 
triggered numerous proposals in which the cosmological 
constant is exchanged with scalar field sourced dynamical 
Dark Energy |7j or alternative theories of gravity 

The set of viable theories is severely limited by Ostro- 
gradski's Theorem [3]. It states that there exists a lin- 
ear instability in any non-degenerate theory whose fun- 
damental dynamical variable appears in the action with 
higher than 2nd order in time derivatives: the Hamilto- 
nian for this type of theory is not bounded from below 
and therefore it accepts configurations with arbitrarily 
large negative energy [TOl [11]. This result can be by- 
passed by considering degenerate theories, i.e. those in 
which the highest derivative term can not be written as a 
function of canonical variables. In this case, the dynam- 
ics is described by second order equations of motion, even 
while the action contains higher derivative terms. If grav- 
ity only involves a rank two tensor. Lovelock's Theorem 
[l2] states that the Einstein-Hilbert action with a Cos- 
mological Constant is the only theory based on a local, ^ 
Lorentz-invariant Lagrangian depending on the metric 



^ There is some evidence that there exist viable nonlocal theories 
ameliorating the fundamental problems of gravity I13| . 



tensor and its derivatives which gives rise to second or- 
der equations of motion in four space-time dimensions. 

The addition of a scalar degree of freedom provides a 
generous extension of the possibilities. The most general 
gravitational sector for a scalar-tensor theory was first 
derived by Horndeski [14] and has received considerable 
attention recently |15H21j . It is given by the Horndeski 
Lagrangian 
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Up to total derivative terms that do not contribute to the 
equations of motion, the different pieces can be written 
as M 
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Here R, G^i, are the Ricci scalar and the Einstein tensor, 
X = —^g^^^(j)_^(j)^f^ is the scalar field canonical kinetic 
term and commas and semi-colon represent partial and 
covariant derivatives respectively. On top of a generalized 
k-essence term ([2]), the remaining pieces (3][5 1 fix the ten- 
sor contractions, which rely on the anti-symmetric struc- 
ture of the terms to trade higher derivatives with the 
Riemann tensor in the equations of motion. Note that 
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Einstein gravity is recovered by a constant G4 = Mp/2, 
while a field dependence G4 = oj{(p)Mp/2 yields an old 
school scalar-tensor theory, without adding higher deriva- 
tive interactions (when combined with a suitable kinetic 
term for the scalar e.g. Brans-Dicke [22j). The theories in 
which the free functions in (3][5) depend on the canonical 
kinetic term X require the presence of degenerate terms 
with higher derivatives. Theories for which G^^GijGt^ 
have simple X— dependences are usually known as covari- 
ant Galileans |15H17j . while theories with more general 
X dependence are often known as generalized Galileans. 
Some of the possibilities considered so far are listed in 
Table H 

In the pursue of generality, one can further consider 
theories in which the scalar field is allowed to enter the 
matter sector directly. This type of relation is found in 
old school scalar-tensor theories, which can be expressed 
as Einstein's theory, plus a scalar field entering the mat- 
ter sector by means of a conformal transformation j8j. 
Bekenstein studied the most general relation between the 
physical and the gravitational geometry (i.e. the two 
metrics out of which the gravitational and the matter 
Lagrangians are constructed) compatible with general co- 
variance [55]. When it only involves a scalar field 0, it is 
given by the disfarmal relatian 



(6) 



The free functions A and B may also depend upon the 
scalar kinetic term X in general, but we will focus on the 
simpler case here. Previous applications of such a rela- 
tion to cosmology include varying speed of light theories 
[2^ [55] , Lorentz invariance violation [3S] , inflation [57] , 
massive gravity [28], Dark Energy [29l [30], relativistic 
MOND theories (311 [32] and extensions of Dark Matter 
[551 The present work studies the implications of 

such a coupling, expanding on our previous analysis [35] . 

The disformal relation ^ can also be motivated in the- 
ories with extra dimensions, in which matter is confined 
to a 3-1-1 dimensional brane embedded in a larger bulk 
space [361146] (see Ref. [47] for a review). The action 
for this type of theories is constructed using geometric 
scalars computed out of the induced metric 



(7) 



where the moduli fields tt^ represent the coordinates or- 
thogonal to the brane and g^i, is the bulk metric prior 
to the embedding, necessary to describe gravity. In the 
case of a single extra dimension |36], the most gen- 
eral Lagrangian contains four terms with a particular 
form of the Horndenski free functions ([2][5| and arbi- 
trary prefactors. The quadratic term is due to the brane 
tension and has the Dirac-Born-Infeld (DBI) [48] form, 
G2 oc \/^ + (97r)2. Therefore, these models are known 
as DBI Galileons [36' . The higher order terms arise from 
curvature invariants computed out to the induced metric 
([7]) , which produce second order equations of motion |12| : 
G3 arises from the trace of the extrinsic brane curvature. 



G4 from the Ricci scalar and G5 from a combination of 
extrinsic curvature terms and the induced Einstein ten- 
sor. DBI Galileons with more than one extra dimension 
only accept the generalization of the quadratic and quar- 
tic terms G2 , G4 in their Lagrangians [571 ESI • This re- 
striction is necessary to preserve the symmetry between 
the directions transverse to the brane - e.g. the moduli 
fields TT/ in Q. 

The usual Galileon terms [49] are obtained from DBI 
Galileon Lagrangian by assuming a flat bulk metric 
.9/ii' ~^ Vfiv ^-iid taking the non-relativistic limit (i.e. low 
order corrections in {Ott)'^). Galileon theories have at- 
tracted considerable attention recently [501173] because 
they capture interesting features of higher dimensional 
models such as DGP [73], including the Vainshtein 
screening mechanism [75]. This effect hides the presence 
of the scalar force due to the non-linear derivative self- 
interactions of the field, which suppress the field's spa- 
tial gradients around matter sources within the so-called 
Vainshtein radius. The extra force is active on larger 
distances, potentially having significant cosmological im- 
plications. 

This work presents results that may simplify consider- 
ably the analysis of theories based on higher dimensional 
models. Section jn] p resents the coupling to point parti- 
cles and Section ITnl introduces theories constructed out 
of two metrics which are disformally related. In Section 
|III A| it is argued that by performing a disformal trans- 
formation, a theory in which the gravitational sector is 
standard, but the matter metric is constructed disfor- 
mally can be rewritten in a form equivalent to the 
quartic DBI Galileon Lagrangian Q, which arises from 
the scalar curvature computed using the induced met- 
ric Q [36 . Disformally coupled theories therefore pro- 
vide an Einstein Frame description of certain brane- world 
constructions, similar to the way in which the field de- 
pendent coefficient of the Ricci scalar can be moved from 
the gravitational to the matter sector in old school scalar- 
tensor theories. In Section |IV[ the equations of motion 
are derived in the Einstein Frame, in which the gravita- 
tional sector has the Einstein-Hilbert form but the matter 
action includes the scalar field as prescribed by Equation 
([6]). Some properties of the field and the coupling are 
discussed in Sections IIV Al and IIV Bl 

The Einstein Frame description of disformally coupled 
theories unravels the existence of a disfarmal screening 
mechanism [35! , in which the coupling vanishes if the 
field is static and the coupled matter behaves as non- 
relativistic dust. Section jVj explores the dynamics in 
high density, non-relativistic environments for the sim- 
pler case of a canonical scalar with a potential and no 



conformal coupling. In Section VA a simple solution is 
derived in which the field rolls homogeneously regardless 
of the matter distribution, hence avoiding the formation 
of spatial gradients that would give rise to an additional 
force. Section VB presents some regimes in which the 
effects of the disformal coupling might be observed. 
Due to the equivalence between the disformally cou- 



3 



Theory 


Ga 


G3 


(j4 






General Relativity 


A 





A/2 
2 







Quintessence 
Generalized k-essence* 


X + V{<f>) 
K{X,cb) 






Afp 
2 

2 








Old school Scalar- Tensor: 

- Jordan Frame 

- Einstein Frame 


X + V{(j}) 

X + v{4>) 






Ml 
2 








Covariant Galileon^ [15] 
Kinetic Gravity Braiding [76] [77] 
Purely Kinetic Gravity [78] 


K{XA) 
X 


-^X 

G{XA) 



2 AfB 

2 

Ml 
2 


3c5 Y"2 




guv 


DBI Galileon iMj 










gt^i^ 


Disformally Coupled Scalar [271 [35] 


X + V(0) 














^ See Table 1 in Ref. |30| for an assortment of k-essence models constructed using disformal relations. 

§ The usual Galileon |49| is recovered in the absence of curvature. The analysis of these theories often postulates a conformal coupling 
between the matter and the field (conformal Galileon). 

t References |41l I42| provide generalizations constructed within the probe brane scheme. 

Table I: Horndenski projection of Modified Gravity and Dark Energy theories. The possibilities shown take into account 
the Horndenski Lagrangian ([T| (see also |19j). the argume nts by Bekenstein leading to the disformal metric ([gJi and the 



possibility of defining disformally related frames (see Section III I. It is then possible to consider a theory of the form Sub = 
/ d'^x (^^—gCn + \f--g^ J^Mig'^vi^)^ as the most general case with a universal coupling to matter. Here Mp = (87rG)~^, 
X = — |0,n<?!>''' and 7 = ^-^^^x '^^ ^ brane Lorentz factor. 



pled theory and a covariant Galileon, the disformal 
screening effect might just be an alternative way to re- 
gard the Vainshtein mechanism and study its properties. 
These screening mechanisms rely on the higher deriva- 
tive form of the field kinetic terms (Vainshtein) and the 
kinetic mixing between the field and the coupled matter 
(disformal). The other available screening mechanisms 
are essentially different, as they exploit the interplay be- 
tween the field potential and the coupling to matter: the 
chameleon fields rely on the high mass of the field in dense 
surroundings [79|, and the symmetrons are screened in 
high ambient density due to their field-dependent cou- 
pling I8OJ . Screening mechanisms are central to the con- 
struction of alternative theories of gravity in which mod- 
ifications are allowed to occur over cosmological scales, 
while the gravitational physics operating in the Solar Sys- 
tem are close enough to GR to satisfy current precision 
tests [HI]. 

The cosmological implications for these models are 



considered in Section VI The intensity of the purely 
disformal coupling is approximately proportional to the 
scalar field energy density p^, unlike in the conformally 
coupled case for which it is proportional to the coupled 



matter density pm- The equations for non-relativistic 
coupled fluids are given at both the background and 
linearly perturbed level, including the analysis of fixed 
points and an analytic expression for effective gravita- 
tional constant on small scales. The cosmological equa- 
tions are solved numerically for a simple Disformally 
Coupled Dark Matter model (DCDM), in which Dark 
Matter is the only coupled species. The model is pre- 
sented and analyzed in detail in Section |VI A[ includ- 
ing the computation of Dark Matter and baryonic power 
spectra. The simple DCDM model enhances the growth 
of the coupled Dark Matter density contrast too much 
to be compatible with observations. However, the model 
contains significant freedom to provide phenomenologi- 
cally successful alternatives. Several possibilities to ren- 
der the model viable are discussed in Section [VlBI 

We conclude in Section IVIII with a discussion of the 
main results and future research directions. Appen- 
dices [5|D] contain a summary of disformal relations and 
some lengthy expressions that were not necessary for the 
main discussion. Throughout the present work, quan- 
tities computed or constructed using the metric ^ are 
denoted "barred" or "disformal" . Quantities constructed 
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using are denoted as "unbarred" and do not involve 
the scalar field. The metric signature is (— , +, +, +) and 
units in which the speed of light c — I are assumed unless 
specified otherwise. 



II. A TEST PARTICLE IN A DISFORMAL 
METRIC 



Let us start with the simple exercise of determining the 
dynamics of a point-like particle with mass m coupled to 
the disformal metric ([6]). A Lagrangian density for such 
a system is given by 



gCp = -m^^g^,x^x-5''^\x^' - a;^(A)) , (8) 



where the dot means derivative with respect to the affine 
parameter A along the trajectory a;(A) and the correct 
weight for the delta function has been taken. ^ The ef- 
fects from the coupling can be seen from the barred four- 
velocity modulus in (Isl) 



(9) 



Distances are dilated by the conformal factor A, as usual. 
The disformal factor B gives an additional direction- 
dependent effect proportional to the projection of the 
four- velocity along the field gradient. The equations of 
motion can be obtained by maximizing the proper time 
of the particle along its path. The result is the disformal 
geodesic equation 



0, 



(10) 



where the barred Levi-Civita connection has been as- 
sumed to be torsion-free and such that the metric 
compatibility relation holds for barred quantities, i.e. 
Vq^/ji/ = 0. It can be computed from (|6| and written 
in terms of unbarred covariant derivatives of the barred 
metric in a rather compact form 
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Here the symmetrization is defined as t(Q/3) = 
\ {tap +ipa)- No assumption about the dependence of 
A, B has been made to obtain the above expression, 
which remains valid if A,B depend on X. Note that 
the difference between the two connections is a tensor. 



as expected. Appendix B shows the expansion of (11) in 



terms of A, B and its derivatives, which is rather lengthy 
to be included here. 

The stress energy tensor with respect to the unbarred 
metric can be computed by variation of (|8| with respect 
to g^i, 



2 5{^Cr. 
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If the gravitational metric is the unbarred one, this is the 
energy momentum tensor sourcing the space-time geom- 
etry.'^ This result can be used to express the particle 
Lagrangian in terms of the energy momentum tensor 



- -g^^T^'' . (13) 

The above expression gives an effective form for the cou- 
pling to matter. It shows how the kinetic term of the 
scalar mixes with the matter content, a very important 
property that lies at the heart of disformally coupled the- 
ories, including the disformal screening mechanism ex- 
plored in Section [V] 



III. DISFORMALLY RELATED THEORIES 

The previous Section presented a simple example of 
a theory in which the matter Lagrangian is constructed 
using a disformal metric ([6]). Although no gravitational 
sector was specified, the simplest possibility is to assume 
that it is given by the Einstein-Hilbert form computed 
out of the unbarred metric 17^^. In this case, Einstein 
equations retain the usual form and are sourced by the 
energy momentum tensor ([l2|. We shall refer to disfor- 
mally coupled theories in which the gravitational sector is 
standard as being expressed in the Einstein Frame (EF) , 
in analogy with old school scalar-tensor theories. More 
generally, one wishes to know what kind of theories can 
be constructed using two metrics that are disformally 
related and study the connections between them. This 
generalizes the conformal equivalence between old school 
scalar-tensor theories minimally coupled to matter and 
theories with a standard gravitational sector, but with a 
non-minimal coupling between matter and the scalar. 

In order to consider theories which allow an Einstein 
Frame description, one starts with a general bi-metric 
theory where the gravity sector has the EH form, but 
with unspecified forms for the gravitational and matter 



^ The one-dimensional definition of the delta function requires that 
its generalization to higher dimensions cancels out the tensor 
density in the integrand S^j^\x — xq) = -y=IlaS£){x°' —Xg) (e.g. 

in spherical coordinates (r^ sind)~^S^^\x) = 5£,{r)5jj{d)5£){(f))). 
Hence it does not matter whether \/— g or y/—g is used in the 
integration, as long as the delta function is consistent with it. 



^ It is possible to write ( |12| in the perfect fiuid form T*^" = 
pu^u" if the coupled matter four velocity and the en- 
ergy density are identified with = x^/\/—x'^ and p = 
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Table II: Physical frames for disformally coupled theories. The intermediate frames are named after the effects of the disformal 
coupling, e.g. a disformal coupling in the matter sector (Disformal Frame) and a quartic Galileon term in the gravity sector 
(Galileon Frame). The transformation rules (solid arrows) are based on the action ( |14[ | and given in terms of the definitions 
( 15|18 1. Note that the transformations commute. Dotted arrows indicate whether the conformal and disformal parts of the 
coupling enter the gravitational or the matter sector for each given frame. 



metrics 



(14) 

Playing with the disformal relations between g'^^^ and g^{, 
allows one to write the above theory in different frames. 
Besides the Einstein Frame, an obvious possibility is to 
consider the Jordan Frame (JF), a description in which 
matter appears minimally coupled and the field only en- 
ters the gravitational sector. But since the disformal 
coupling has two parts, two more intermediate frames 
can be defined, in which only a certain part of the cou- 
pling enters the matter action. The four possibilities are 
described below and summarized in Table [IT] together 
with the transformations that provide the connections 
between them. For the sake of simplicity, the Einstein 
Frame has been defined using a matter metric of the form 
([6]), consistently with the notation used in most of the 
paper. 



1. Einstein Frame: 



= 9tiu, g^i = Ag^^ + Bcjy^^cj)^^ . (15) 

This is the formulation used throughout the rest of 
the paper. The equations in this frame are derived 
in Section [ll] for a point particle and in Section |TV] 
in general. 



2. Disformal Frame: 

1 
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9t^i^ + 



(16) 



The disformal part enters the matter Lagrangian 
explicitly. The conformal factor enters the gravita- 
tional sector through a coupling to R, like in old 
school scalar-tensor theories. 



3. Galileon Frame: 

B 

a' 



9t^u 



Agf,^ . 



(17) 



The conformal part enters the matter Lagrangian 
explicitly and the field couples directly to gravity 
as a DBI Galileion, see Section IlII A 



4. Jordan Frame: 



G _ I B 
9i_iu ~ ^5^11^ 
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M _ 
9av ^ 9fJ.iy ■ 



(18) 



Matter is minimally coupled to a metric and the 
field enters the gravitational sector exclusively. 

The JF is the most convenient frame to analyze certain 
properties of the theory and its predictions, as matter fol- 
lows the geodesies of the simple metric g^j^. The matter 
metric in the remaining frames contains the scalar field 
explicitly, and therefore matter moves along geodesies 
that involve the field variations ( 10 1 in these representa- 
tions. These frames are still interesting to analyze the 
theory. For example, the equations simplify considerably 
in the EF, just like in conformally related theories. Once 
these are solved, the solutions can be used to write down 
the Jordan Frame metric. 

The explicit computation of the curvature scalar for 
a metric which includes a disformal part allows one to 
connect the theory studied in the Einstein Frame with a 
particular sector of the Horndenski Lagrangian 0. As 
anticipated in the introduction, the thus obtained the- 
ory is related to a type of DBI covariant Galileon when 
expressed in the Galileon or Jordan Frames. 



A. Disformal Curvature: The Galileon Frame 

It is possible to get a sense of disformally coupled the- 
ories in a different frame by applying the transformations 
sketched in Table [TT| to known actions. The simplest case 
is the canonical scalar field as was described in Ref. [3D] , 
where it was shown that the transformation produced the 
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disformal quintessence Lagrangian 



-g{X ~ V) 
X 



(19) 



, - VA - 2BXV ] , 

^/A - 2BX J 



up to the ambiguity in the definition of the kinetic term 
described in Section 2.4 of Ref. [3^^. The above theory 
encompasses several Dark Energy models in certain lim- 
its, which are obtained by appropriate choices of A, B 
and V. 

Considering similar relations when the disformal trans- 



formations involve the gravitational sector in ( 14 ) re- 



quires the computation of the Ricci curvature for a barred 
metric that includes the scalar field as in ([6]). The start- 
ing point is the difference between the standard and the 
barred connection (11) 



p. ^ - 5"^ (V(^.g.)A - ^ Va5m.) (20) 



where the symmetrization is defined as t(ap) = 
\ {tafi -\-t0a). The barred Riemann tensor is obtained 
from the usual definition, and it can be related to the 
unbarred one in a manifestly tensorial form in terms of 



(21) 



where anti-symmetrization is defined without the usual 



I coefficient = 
from the contraction 



A, 



- Abo ■ The Ricci scalar follows 



with the inverse barred metric 



^ AV A-2BX^ ^ 



(22) 



(23) 



Finally, the disformal Einstein-Hilbert Lagrangian den- 
sity requires the barred volume element (A4| to be co- 
variant 



(24) 



Note that no assumption has been made about the func- 
tions A, B out of which the geometric quantities ( 20p2 ) 
are computed. However, the general computation is very 
lengthy, and it is useful to adopt some simplifications. 

Let us focus for the time being on a theory in the 
Galileon Frame, for which the disformal part is absorbed 
into the gravitational sector. The following computation 
assumes thus a gravitational metric of the form (17) 



9^ 



D 



(25) 



where the disformal factor D{(j)) is yet to be specified.^ 
The Jordan Frame can be obtained at the end of the 
computation by inverting the conformal transformation 
5^1/ — >■ A~^((l))gfj_^ in the resulting metric and curvature 
objects. Since the transformation rules for curvature ten- 
sors under conformal relations are well known [82], and 
Galileon-like theories usually retain a conformal coupling 
to matter in phenomenological applications (e.g. 83 ), 
the Jordan Frame curvature will not be computed ex- 
plicitly. 

Let us further redefine the field to simplify the barred 
metric (25) 



The above expression has the same form as the effective 
metric in probe-brane theories ([t]). It simplifies consid- 
erably the computation of the connection tensor ( 20 1 , 
which now reads 



Here 



1 



7 = 



(27) 



(28) 



^l + vr^^TT'" ' 

is a Lorentz factor that arises from the inverse met 



ric (23). The barred Riemann tensor can be com- 



puted directly from (21), rewriting anti-symmetrized 
derivatives in terms of the curvature V[^Viy]X'.'.'^'.'.'.' = 

- -Sji?^^^.^^^ :.^:..- The resuh is also sim- 
ple 



(29) 



The barred Ricci scalar can be easily obtained by a sec- 
ond contraction 

R = {g^, - 27'7r,^7r^,) [^"^+72 (tt-^-^D^ - ^''^^tt'"'^) 

(30) 

Finally, the gravitational Lagrangian y/^gR ^ust requires 
multiplying by the barred volume factor ( 24 ) to make it 
covariant. 

The total action for the theory in the Galileon Frame 
is obtained by adding a matter Lagrangian with a con- 
formal factor A{(j)) in the matter metric is 



-g£„i) = / d'^Xy/^C 



GF 



(31) 



* Just as in the rest of the paper, D has been assumed to be inde- 
pendent of the field derivatives X. This assumption is important 
in order to simpUfy the computations performed below. However, 
given the importance of the X dependence of higher derivative 
terms in the Horndenski Lagrangian ([TJ and related theories, it 
is worth considering the more general case D{rp,X), although in 
this case the equations would become very involved. 
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with 



GF 
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(32) 



These results were previously obtained by de Rham and 
Tolley |36| in the context of higher dimensional grav- 
ity theories. The above action corresponds to a quartic 
DBI Galileon term, with the right form of the Horndcn- 
ski Lagrangian Q after integrating by parts. It reduces 
to the quartic covariant Galileon term in the non rela- 
tivistic limit limit X = — ^Tr^^TT'^ <C 1, with the right 
non-minimal coupling to gravity. If the number of extra 
dimensions is larger than one, only the above and the 
quadratic DBI brane tension terms are allowed, as they 
respect the symmetry between the directions transverse 
to the brane [37l |43j . The reader is referred to Section 
5 of Ref. [3 6) for results on other curvature invariants 
leading to the cubic and quintic DBI Galileon terms for 
probe branes with a single extra dimension. 

The Jordan Frame expression of covariant Galileons 
was considered by Appleby & Linder [SS], where both 
the conformal and the disformal coupling to matter were 
shifted to the gravitational side of the action. The the- 
ory studied there is essentially different from ours, as 
it includes all the standard Galileon terms (which are 
obtained from the DBI Galileon terms in the limit of 
small X), instead of considering only the quartic term 
([32| ). In the context of inflation, Renaux-Petel et al. 
[oSl [55] studied a theory including quartic DBI Galileon 
terms in the Jordan Frame. Their model also included 
an Einstein-Hilbert term computed out of the unbarred 
metric. Therefore, it does not allow for the construction 
of an Einstein Frame and essentially differs from the the- 
ory considered here. However, their study stresses the 
bi-metric structure of such theories and uses interesting 
techniques to analyze the dynamical equations. 



IV. EQUATIONS IN THE EINSTEIN FRAME 



It will be further assumed that there is no dependence 
on the barred metric derivatives.^ A scalar field La- 
grangian density of the k-essence type — C^{(j), X) has 
been also included. More general dependence on the field 
derivatives may be considered, but this term gives rela- 
tively simple equations of motion, and is general enough 
to accommodate both a canonical = X — V and a dis- 
formally self- interacting scalar field [221 ES] ■ The matter 
Lagrangian may include other pieces with different cou- 
plings. An uncoupled matter sector can be included by 
the addition of a Lagrangian \J~gl^u constructed out of 
the unbarred metric. 

The stress energy tensors for both species will be fur- 
ther defined in terms of the contravariant gravitational 
metric 



rji(J,V 




^ Sg, 



(35) 
(36) 



such that the Einstein field equations take the usual form 
G^" = SttGT'^'^ and the total energy-momentum is co- 
variantly conserved with respect to the unbarred metric 
by virtue of the Bianchi identities: V^(T;^'' -I- T^") = 0. 
However, the coupling causes that this relation does not 
occur for each component separately, and in general 



(37) 



where the form of the interaction (last equality) can be 
seen by explicitly computing the divergence of the scalar 
field stress tensor: 



= c 



(38) 



As the term in parenthesis is equal to the Lagrangian 
field variation 5C^/54>, the equation for the scalar field, 
SS/S(j) = 5C^/5(j) + SCm/6(j) = 0, allows one to write Q 
in terms of the variation of the matter Lagrangian 



Q 



(39) 



In this Section, the equations for a disformally coupled 
theory which admits an Einstein Frame description will 
be derived. Such a theory is given by the following action 



'EF 



I&ttG 



-g^mig^L,^,ip) 



(33) 

The interacting matter sector •J—gCm is to be con- 
structed using the barred metric ([6]), as it has been made 
explicit in the second line using the explicit form of g^v 
and its determinant ( 24 1 



-gt^m{gtiu,i') = 



= ^-gA\ll- 2^A£™(Agp. 



(34) 



The coupling can be evaluated by application of the chain 
rule. For the specific form of the barred metric ([6|, the 
variation with respect to the field yields 



^ iV^^m) dg^,^ dg^j,^ _ ^-g 
hi^y dg^Li' d(f> 2A 



Ti^' (A'5^, -f B'b^K) 



(40) 



This implies that the disformal connection does not appear 
in the action. This assumption holds for scalar fields and gauge 
vectors vanishes due to the lack of indices and antisymmetry of 
the kinetic term, respectively. Although other fields may couple 
to g^v^x, the assumption simplifies the equations considerably 
and is common in the analysis of scalar-tensor theories. 
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and similarly for the derivative with respect to the field 
gradient 




Sgp.i' dg^,^ d(t) 



The coupling can be obtained after replacing ( 40|41 ) 
in pgl), 



2A 



(42) 



The equations for the coupled matter component and the 
field are then 



2XC 



4>,X4> 



Q, 



(43) 
(44) 



where A^J'^) 



iCp.xg'"' + C^.xx(l)'^<P''') is the general 
kinetic term for the scalar. Einstein equations G^" = 
87rG(r^'' + r^'^^p together with (|43| [il and (III) deter- 
mine unambiguously the evolution oimatter, tne scalar 
field and the metric. These equations naturally contain 
the case of a conformally coupled field, where only the 
coupling to the trace of energy momentum is present in 
Q. Note that so far this result is general and does not 
depend upon the matter content as long as the matter ac- 
tion only depends on the field through the barred metric 
^ algebraically. 



A. Properties of the Scalar Field Equation 



The first term in the coupling ( 42 1 contains higher 



derivatives of the variables T^'^.^,^^^ due to the kinetic 
mixing in the mater action, cf. (13). These have to be 



solved for in order to integrate the evolution equations 
(43 44), which can be done after adopting a coordinate 



system. It is possible to eliminate the matter derivatives 



in the scalar field equation by contracting ( 43 ) with 
and solving for 

back in (44) and rearranged as 



y ^Tif^ . The result can be inserted 



■M^-V^V.^ + ^ Q^.rr + V ^ , (45) 

were we have defined: 

M>^- = C^,xg^^ + C^^xx^^^r - , (46) 



A-2BX 



_ A! (MB B' 

V = £ rf, - 2XC 



:X4 



(47) 
(48) 



This equation can then be used instead of (44) to de- 
termine the evolution of the scalar field. It displays very 
clearly the role of the coupling, which enters not only as a 
modification to the effective potential (second term) , but 
also in the coefficient for the higher derivatives of the 



field. This feature will be ultimately responsible for the 
screening mechanism that these models exhibit in high 
density regions, which is explored in Section [V) 

Equation ( 45 ) also shows that the different components 



of the energy momentum tensor may modify the signs 
of the coefficients of the second order field derivatives. 
It is then necessary to analyze whether the scalar field 
propagation has a good initial value formulation, as was 
pointed out and first ana lyze d by Bruneton & Esposito- 
Farese [84].^ Equation (45) is a quasi-linear, diagonal 
second order equation [8o|. of the form 



(49) 



However, its hyperbolic character relies on the signa- 
ture of A^'^", which involves the coupled matter energy- 
momentum tensor. For a canonical scalar field = 
X ~V disformally coupled to a perfect fluid, the deriva- 
tives tensor reads M'^^ = a-2BX '^^^si~PTPjPTP) in 
coordinates comoving with the fluid. Positive energy den- 
sity keeps the correct sign of the time derivative term if 
B > 0, avoiding the existence of ghosts modes. However, 
a large pressure can flip the sign of the spatial derivatives 
coefficient, introducing a gradient instability. This might 
have important consequences in sufficiently relativistic 
environments. 

Addressing the viability of the theory hence requires 
determining under which conditions the instability may 
occur dynamically, which in turn requires considering the 
evolution of the coupled matter components including the 
non- linear terms in ( 42 45| . In certain cases, the system 
might respond to a situation in which Bp ~ 1 — BX by 
diluting the (Einstein Frame) pressure below the thresh- 
old value or softening the spatial gradients of the scalar 
field. In this sense, the instability induced by the pres- 
sure may be analogous to the potential existence of sin- 



gularities in the disformal volume element ( 24 1 whenever 
g (X A — 2BX — > 0. This singularity in the barred metric 
is avoided by the field evolution, as it slows down when- 
ever B4>^ — > A. The mechanism exploited to induce a 
slow roll phase in cosmological applications is precisely 
this dynamical resistance to pathology (cf. disformal cou- 
pling to matter described in Section fVI A| and disformal 
quintessence p9l l30] ). 

Studying the conditions under which the pressure in- 
stability can be avoided dynamically might as well re- 
strict the allowed functional forms of the conformal and 



® Disformally coupled theories are discussed in Section IIIC of ref. 
84, from a field theoretical perspective. There it is argued that 
some dependence on X on either the conformal or the disformal 
factor is necessary for the theory to have a Lorentzian signature 
for all values of X (cf. their eq. 3.22). However, the dynam- 
ics of the scalar field might prevent it from acquiring arbitrarily 
values of X, as it is argued below and shown explicitly for a cos- 
mological model in which 2BX < 1 at any time (See Figure 1). 
Ref. 1841 also studies the hyperbolic condition within dynami- 
cal pressureless matter, which is complementary to the analysis 
presented below. 
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disformal factors. In the worst case, it might spoil the 
disformal screening mechanism, or even completely for- 
bid the occurrence of a disformal coupling. Determin- 
ing whether or not this is the case will be the objec- 
tive of future work. The Einstein Frame pressure p will 
be neglected as a part of the approximation scheme in 
the following analysis, implicitly assuming that Bp <C 
A — 2BX, in order to avoid potential issues. 



B. Coupling to Perfect Fluids 



Assuming a perfect fluid in the Einstein Frame T^'^ — 
[p + p)u^u'^ + pg^"^ with u^Ua — —1, it is instructive to 



project (43 1 along and perpendicular to the matter four 
velocity. This determines how the local law of energy 
conservation and the geodesic equation are modified by 
the coupling 



(50) 



for static, pressureless configurations [57]. This is obvi- 
ous from ( 43 ) , since only the component is nonzero 



for dust, and when contracting the field derivatives with 
the stress tensor, a non-vanishing result requires time 
evolution of the scalar field. Therefore, addressing the 
effects of disformal couplings requires studying the field 
dynamics in high density environments.^ As it will be 
shown in the next subsection, the kinetic mixing induced 
by the disformal coupling makes the scalar field evolu- 
tion insensitive to the matter distribution, as long as it 
is sufficiently non relativistic and its energy density is 
high. 

For the sake of concreteness, let us restrict ourselves 
to a canonical scalar field coupled to a perfect fluid. The 



general equation ( 45 ) then reads 



(52) 



fA' 


fA'B 








[ A 


t) 





(p+p)w"V„u^H-[5^" (V„p+QV„0) = 0. (51) 

In the first equation the coupling modifies the energy con- 
servation relation, due to the energy transfered from the 
scalar field, which is modulated by the projection of the 
field gradient along the 4-velocity. The second equation 
determines the departure of geodesic motion with respect 
to the gravitational metric. The first term describes the 
force arising from the pressure gradient and the second 
the additional force exerted by the scalar field. Both 
forces are projected into the direction _L (coefficient 
in brackets) due to the orthogonality of the four velocity 
and four acceleration. 



The analogue of ( 45 ) for the covariant matter conserva- 
tion equation without second order field derivatives can 
not be obtained without choosing a time slicing due to 
the different high derivative structure in both equations. 
Nevertheless, there is no need to do so, since we already 
found a bona fide field equation ( 45 1 that can be inte- 



grated consistently with the corresponding equation for 
matter (43), substituting the appropriate value of Q. It 



is possible to solve for the time derivatives of all the vari- 
ables after a metric ansatz has been chosen, as will be 
done in Section |Vl] for the study of FRW models and 
cosmological perturbations. 



where X = A — 2BX < A is bounded in order to avoid 
a singularity in the vol ume element of the barred metric 
cx A"^/^, cf. (24). The form of the field equation 



strongly suggests that the dynamics of the coupled sys- 
tem will be different in high than in low density environ- 
ments. Since the energy momentum-tensor appears as a 
coefficient of the higher derivatives as well as in the effec- 
tive potential, there is a well defined limit T*"^ = p oo, 
in which the field equation simplifies considerably. This 
property will be crucial for the disformal screening mech- 
anism 1351. 



A. Dense, Non-relativistic and Static Matter 

The study of Solar System and laboratory tests of grav- 
ity requires considering energy densities that are much 
higher than the cosmological average and pressure is com- 
pletely sub-dominant. As a first approximation, this 
regime can be explored using the general scalar field 
equation (52) for a static, non-relativistic matter distri- 
bution p{x) in the limit p 



oo. More precisely, the 
following dimensionless ratios will be assumed to be neg- 



V. THE DISFORMAL SCREENING 
MECHANISM 

In this Section we consider the compatibility of grav- 
itational theories based on a disformal coupling and lo- 
cal gravity tests. Due to the stringent bounds on scalar 
forces and post-Newtonian effects [SH [86] , some sort of 
screening mechanism is necessary to hide the coupling in 
dense environments such as the Solar System. The disfor- 
mal contribution to the conservation equations vanishes 



''' It is possible to obtain some insight into the dynamics of the field 
from the analysis of the background cosmology given in Section 
|VI| where it is argued that for the purely disformal case the cou- 
pling was proportional to the scalar field energy density, cf. Eq. 
( |63| . This causes the existence of two regimes, a matter domi- 
nated regime in which the effects of the coupling are small, and 
a field dominated regime in which the coupled matter equation 
of state is modified, see Figure ^ When denser regions form, 
the scalar field energy density becomes insufficient to produce 
large effects on the matter distribution, unless the field gradi- 
ents follow the matter distribution and intensify the additional 
force. 
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ligible 




X X 

' "Bp' lip 



V'/'^, r^o-^.p/f^ ^0, (53) 



and Bp < A~2BX, as argued in Section [lV A[ These ap- 
proximations quantify to what extent the effects of grav- 
ity and pressure are disregarded, the requirement of hav- 
ing "soft" spatial gradients relative to the time evolution 
and the fact that Bp is large enough. These conditions 
will be briefly discussed at the end of the Section, focus- 
ing on systems in which their lack of fulfillment might 
lead to observable signatures. 



tion (52) considerably 



The set of assumptions ( 53 ) simplifies the field equa- 




2M„ 



where the first equality is general and the second applies 
to a purely disformal coupling with exponential forms, 
such as the example cosmological model presented in 
Section [VI A[ The above expression departs substantially 
from the simple conformal coupling, for which the p — > oo 
limit is ill-defined. Two important features of the above 
equation endow the theory with the disformal screening 
mechanism: 

• The spatial derivatives become irrelevant, as they 
are suppressed by a, p/p factor with respect to the 
time derivatives.^ 

• The equation becomes independent of the local en- 
ergy density, making the field evolution insensitive 
to the presence and distribution of massive bodies. 



proportional to (j), stronger couplings decay earlier. It is 
possible to obtain a solution for A = 1, p = keeping the 
potential V, but otherwise assuming the simplifications 
(53 1. It is given as an implicit function 



t-to 




2V{c^')/p 



(56) 



where (i^ = % 



^ — ^^^-^ The potential appears sup- 
pressed with respect to the energy density. In tracking 
Dark Energy models, such as the one explored in Section 



VIA is a decreasing function of the field and > 0. 
Therefore, V is of the order of magnitude of the average 
cosmic density po £^nd can be safely neglected if p is much 



higher, recovering the simpler solution (55). As the field 



slows down with time if B' /B > 0, the order of mag- 
nitude of the field time derivative is also cosmological, 
(54) <P^Ho. 

One of the effects of the coupling is to modify the 
energy conservation equation for matter in the Einstein 
Frame, cf. Eq. (50), inducing a variation of the gravita- 



tional mass. In a gravitationally bound two body system, 
this effect is degenerate with a possible time evolution of 
Newton's constant § ^ ^'^ ^ ^^'^^ approximation, 

as can be argued by deriving the expression for the New- 
toninan force with respect to time. Lunar laser ranging 
measurements place precise bounds on this effect to the 
level of G/G < 10"^/Gy [SSj. The magnitude of energy 
density variation induced by a disformal coupling can be 
estimated as p ~ — (0 -I- V')4>- Assuming (j)"^ ^ V ^ po a.s 
discussed above, BpQ > 1^'^ and the solution (55), typical 
mass variation rates M/M are as small as ~ 10~^/Gy for 
the interstellar medium and ~ 10~^^/Gy for the average 
Earth density, well beyond the sensitivity of Lunar laser 
ranging measurements. 



These features ensure that the field rolls homogeneously 
and avoids the formation of spatial gradients between 
separate objects, which would give rise to the scalar force. 
The above properties are caused by the kinetic mixing 
between the field and matter degrees of freedom, and lay 
at the core of the decoupling between both components. 

Let us analyze the simpler, purely disformal exponen- 
tial case. The second equality of equation (54) can be 
easily integrated 



(55) 



In this solution the field time variation is approximately 
constant while t <C 



/30(O) 



and slows down afterwards as 



cx 1/t. Since the coupling to non-relativistic matter is 



Equation ([54| also follows from taking the limit p S> A/B, ij)^ in 
the FRW coupling density Qo l |60| l, precisely due to the absence 
of spatial derivatives. 



B. Potential Signatures 

New local, astrophysical and cosmological signatures 
may be found by relaxing the approximations ( 53 ) as- 
sumed in the previous analysis. Some situations where 
the coupling might become observable include: 

a. Matter velocity flows: The spatial component of 
the matter four-velocity T'^^ mixes the time and space 
derivatives of the field, which may source the field evolu- 
tion (as non-zero velocities introduce terms proporional 



Under these assumptions, Eq. 



V' 
Bp 



V2- 



52|l can be written as <f>-\- ^'p'^ + 
0. The second equality can be 



Co 
B 



Bp 



which 



directly integrated, giving the constraint < 
can be integrated again to obtain (|56|. 

This assumptions are made for order of magnitude estimates. 
For cosmological applications in which the field's potential ac- 
celerates the Universe, p, V are significantly larger than 0^ and 
Bp > 1 occurs while B4>'^ < 1, cf. Figure 2 
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to (j)-oi and (t>(t>,i)- These eflFects are suppressed by a rela- 
tivistic v/c factor, but they may be important in certain 
systems such as binary pulsars. 

b. Pressure: Apphcations of the disformal coupHng 
in the context of Dark Energy arguably require a value 
Bpo 3> 1, where po is the average cosmic density. Then, 
even though the pressure is usually negligible with re- 
spect to the energy density, it should be easy to find 
systems for which Bp is also much larger than one. This 
might have important consequences for the stability of 



the theory, as was briefly discussed in Section IV A 



c. Radiation: Unlike in the conformal case, the 
disformal coupling has non-trivial effects on ultra- 
relativistic fields for which T 0, cf. (42 1. Some au- 



thors have initiated the study of the disformal coupling 
in scenarios featuring radiation. Brax et al. |89| con- 
sidered high-precision, low-energy photon experiments, 
which might be able to detect the influence of a disfor- 
mal coupling on top of a conformal one. The distortions 
in the baryon-photon chemical potential induced by a dis- 
formal coupling and their signatures on the CMB small 
scale spectrum have been studied by van de Bruck and 
Sculthorpe [90] . Other effects may follow if Electromag- 
netism is formulated in terms of the barred metric, such 
as varying speed of light or modified gravitational light 
deflection [51]. 

d. Strong gravitational fields: The connection coeffi- 
cient rQQ0 ^ in the field derivative term is not suppressed 
by Bp. It represents the effects of gravity, and was ne- 
glected because it is small in most Solar System appli- 
cations, since Fqq = ^^(r — 2GM) in the Schwarzschild 
metric. However, this term might become relevant in 
strong gravitational fields, such as the vicinity of black 
holes or compact objects. 

e. Spatial Field Gradients: In the Bp ^ 1, p ^ p 
limit, the equation for the scalar field (54) becomes inde- 



pendent of the matter content and the field derivatives. 
Therefore, if the field acquires a spatial modulation be- 
fore reaching this limit, it will be preserved by the sub- 
sequent evolution. Spatial gradients of the field formed 
when the linear perturbation theory is valid would then 
be present today, with their actual value depending on 
the details of the transition between the perturbative 
(e.g. the small scale hmit ( [72| ) and the screened regimes. 
Gradients of cosmological origin might be seen as pre- 
ferred direction effects pointing towards cosmic struc- 
tures when analyzed in the Solar System. Spatial deriva- 
tives of the field may also be important if the field is 
rolling sufficiently slow as to overcome the p/ p factor in 
(151). 



These and other settings might lead to characteristic 
signatures and new bounds for disformally coupled theo- 
ries, which will be investigated in the future. It should be 
also possible to obtain the coefficients of the Parameter- 
ized Post Newtonian approximation, which would allow 
a more systematic comparison to local gravity tests. 



VI. COSMOLOGY 

Having addressed the viability of the theory in the 
Solar System, let us consider its cosmological implica- 
tions. Using the Einstein Frame description, the Fried- 
mann equations have the usual form 



SttG 
AttG 



(p+y+^), (57) 
{p + 24>^ - 2V) , (58) 



but the conservation equations for matter and the scalar 



field have to be computed from (43 45 ) 



p + 3Hp = Qo0 , 

4> + 3H<j) + V' = -Qo , 



(59) 
(60) 



were p is the energy density of the coupled matter com- 
ponent and the background coupling factor reads 



Qo 



A' - 2B{3H4> + y + ^(/)2) + 
2(^A + Bip-4>^)) 



(61) 



after solving away the higher derivatives. In the following 
we restrict to flat space, K — 0. 

At this stage it is possible to understand the difference 
between the pure conformal (B — 0) and disformal {A = 
1) cases by writing (61 ) in terms of the equation of state 



and the scalar field energy density: 



fB' , 



2A^ 



(62) 



2y' 



((1 + W^)ptotP0)^^^ , 



(63) 



where in the pure disformal case it has been assumed 
that Bp ^ 1 ^ B(l>^. This approximation is satisfied 
by the model presented in the next subsection when the 
coupling is active, see Figure [2] The last term in (63) 
represents the contribution from the Hubble term, which 
is subdominant when the slopes of B, V are large. The 
above expressions imply that the conformal and disfor- 
mal coupling between Dark Energy and Dark Matter are 
related to essentially different phenomenological param- 
eterizations, where the interaction is either proportional 
to p = Pdm [92H94] or p^ [95l|96]. 

The equations governing cosmological perturbations 
can be obtained from Eq. ( 39 ) , which can be used to read 



both the disformal matter non- conservation and the field 
dynamical equation. Working in the Newtonian gauge 



(1 + 2$)dt2 -I- 0^(1 - 2*)dx2 



(64) 



avoids potential misinterpretations when swapping be- 
tween different frames, at least when these are related by 
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conformal transformations f97' . Solving for the higher or- 
der derivatives, the perturbed continuity and Euler equa- 
tions for the disformally coupled matter contrast S^c — 



5p/p and the divergence of its velocity = ikjT^^.^-^ = 
ikjV^ pa^^ and the scalar field read 



Sdc + - + 



: 



Ode 



= 3^ 



S(j)- 



(65) 



9 H 



(66) 



and 



%^ + {2BV"-B" 
M 



-5Q - 2$ (Qo + V) -f 0(4 + 3*) , 



i-2\ P 



(67) 



2M ' y^B^iy' + ^m + B-^p-^^,^ ^^^^ 



B'(P - B{'SH - pB'(f>) - B^{2V'i + 3i7(p + ^'^)) ^ ^ 



M2 



-B'^ + B{6H - pB'^) + 2B^{3Hp + V'q 



A'P M 



(68) 



r 



for a purely disformal coupling A — 1, where M = 
1 -|- Bp — BcjP' . The coupling perturbation 5Q is given 
in Appendix [C] for the general case. This expression is 
a much more cumbersome combination of the fluid and 
field perturbations than for the purely conformally cou- 
pled case 



^\og{A)'p5d, + \\og{A)"p5cp. 



(69) 



Note that, unlike in the conformal case, the first term in 
(681 is proportional to and hence the coupling intro- 



duces explicit scale dependent terms at the level of the 
equations. This feature will be reflected in the growth 
of perturbations and the power spectrum, studied below 
for an example model. 

To extract the most relevant new features by analytic 
means, we shall consider the subhorizon approximation. 
In the small scale limit, taking into account only the mat- 
ter perturbations and the gradients of the field, there is 
a simple expression for the perturbed interaction SQ. In 
this Newtonian limit, we further relate the field gradient 
to the matter perturbation through the field equation 
(67), which yields the simple expression 



5Q(N) = QoSdc 



(70) 



Combining equations (|65| and (661 together with the 



usual Poisson equation, we obtain the evolution of the 
coupled Dark Matter overdensity 



2H 



Q. 



; 



(5dc = 47rGeffp(5dc 



(71) 



In addition to an extra friction term, the source term is 
modulated. The last effect is captured by defining an 
effective gravitational constant Gcff that determines the 
clustering of Dark Matter particles on subhorizon scales 



Gcs 
~G~ 



Ql 



(72) 



This approximation has the same expression as the sim- 
ple conformal case, although with a significantly different 
functional form of the coupling Qg, which is now given 
by Eq. (1611. 



A. Disformally Coupled Dark Matter 

In what follows it will be assumed that the field is only 
coupled to Dark Matter, while radiation and baryons fol- 
low geodesies of the gravitational metric and do not feel 
the scalar interaction directly - If baryons are also cou- 
pled, then the ratio pdm/pb remains fixed, because both 
species feel the same effective metric. Postulating that 
the baryonic and electromagnetic sectors are constructed 



This can also be seen directly from l |6l| : in the denominator of 
Qo, the energy density has to be substituted by the total one 
p — > pdm + Pb, while the multiplicative coefficient p would refer 
to each individual species. 
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out of the gravitational metric also avoids problems with 
precision gravity tests and the subtleties related to the 
existence of different frames, hence simplifying the anal- 
ysis of cosmological observations. 

To study the dynamics within a particular example, 
we focus on a simple Disformally Coupled Dark Matter 
(DCDM) model, constructed with the following prescrip- 
tions: 

• Dark Matter disformally coupled to a canonical 
scalar field, following Eq. ( 59]|6T |. 



• An exponential parametrization for the disformal 
relation and the scalar field potential: 



V = VQe'^"^'^'" , 
A = I. 



(73) 
(74) 
(75) 



with Mp = (87rG)~i/2^ ^he conformal factor A 
has been set to the trivial value in order to focus 
on the novel features. Furthermore, the coupling is 
chosen to be negligible in the early universe, and 
hence initial conditions and early evolution are not 
affected. 

• Uncoupled baryons, photons and neutrinos, which 
follow the usual barotropic scaling relations p ~ 



a' 



3(l+u)) 



Zero cosmological constant. 



Besides being motivated from some high energy scenar- 
ios, the exponential forms ( 73||75 1 facilitate the choice of 
natural scales for the constant prefactors by shifting the 
zero point of the field (e.g. Bq ~ M~^, Vb ~ M^, Aq 
dimensionless) . Furthermore, these forms allow a conve- 
nient exploration of the phase space of the system. In 
addition to the previously studied fixed points [71 [33] ^ we 
find only one new, a disformal scaling solution that is not 
an attractor. The details of this analysis can be found in 
Appendix [D] 



1. Background Evolution 

The model set up is similar to the uncoupled self- 
interacting field case described in Ref. [3^. In partic- 
ular, the potential ensures a tracking stage for the field 
and the value of <f>Q is chosen to tune the transition time 
when the disformal coupling becomes relevant. Although 
only Dark Matter is affected by the coupling, radiation 
and baryons are included in order to provide a more real- 
istic description. The evolution at early times is then as 
in the usual exponential quintessence model, where the 
field tracks the dominant fluid component and the slope 
of the potential 7 determines the amount of Early Dark 
Energy (EDE) [7] 
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Figure 1: Background evolution of disformally coupled mat- 
ter. Left: evolution of the energy density for the field (red, 
light) and coupled matter (blue, dark) for different choices of 
the coupling slope /?. Right: equation of state for the field 
(red, light) and coupled matter (blue, dark). High values of 
/3/7 (solid, dashed) give a good fit to observations, while low 
values (dotted) do not produce enough acceleration. 



which depends on the dominant matter component equa- 
tion of state parameter Wm- The new features ap- 
pear when the disformal factor BcfP' becomes of or- 
der one. Then the clocks that tick for Dark Matter, 
.900 = — 1 + slow down and make the effective equa- 
tion of state approach minus unity asymptotically. The 
field also slows down to avoid a singularity in the effec- 
tive metric g^^, and the universe enters into a de Sitter 
stage. This natural resistance to pathology was also ob- 
served in the disformal self-coupling scenario described 
in Refs. [29l [30] . The disformal coupling provides then a 
mechanism that triggers the transition to an accelerated 
expansion. The relatively steeper the slope of the disfor- 
mal function is, i.e. the higher the ratio /3/7, the faster 
the transition happens, as seen in Figures[l][2j This tran- 
sition also produces a short "bump" in the equation of 
state, which affects the growth of structure. 

The evolution of Goff for the disformally coupled Dark 
Matter example model ( 73][75 ) is shown in Figure [3j It is 
characterized by a bump at the transition, whose height 
increases with and a further increase when the po- 
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Figure 2: Evolution of the dimensionless disformal factors 
B4>^ (solid), Bp (dashed) and BV (dot-dashed). Higher val- 
ues of /3/a produce a sharper transition and lead to higher B 
at late times. 



tential becomes dominant. At the later stage, the depen- 
dence is approximately Gcs/G—1 ~ {jV/ p)'^ and yields a 
large value since Dark Energy domination requires V > p 
and 7 > 15 is necessary to avoid the effects of early Dark 
Energy ( 76 1 [M] ■ This enhancement occurs on observable 



scales and spoils the formation of large scale structure in 
this particular case. Problematic growth enhancement 
also occurs in conformally coupled models that attempt 
to address the coincidence problem The observ- 

able effects will be analyzed using the full perturbation 
5Q within the disformally coupled Dark Matter model. 
Several alternatives to render the model viable will be 
described in Section IVIB I 



2. Perturbations 



The full system of linearized equations (67]68| was 
solved numerically using a modified version of the Boltz- 
mann code CMBeasy adapted to the Disformally Cou- 



pled Dark Matter model described in Section VI A Since 
matter is essentially uncoupled until z < 10 there was no 
need to modify the initial conditions, which have been as- 
sumed adiabatic. Figure[4]shows the evolution of the den- 
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Figure 3: Effective gravitational constant on small scales ( 72 1 
for different values of /3, 7. The value is large at the transition 
due to the disformal friction term B'ifP' , and latter due to the 
contribution of the potential term BV' (see text and compare 
to Figure [1]). 



sity contrast of disformally coupled matter. The baryons, 
which are uncoupled in this particular example, are also 
shown for comparison. Figure [s] displays the power spec- 
tra for disformally coupled matter and baryons at z = 
for different values of the parameters. Figure [6] shows the 
CMB power spectrum and the baryon-DM bias induced 
by the coupling at z = 0. 

Besides the effect of early Dark Energy and late time 
scalar force captured in Geg, the disformal coupling 
causes a considerable integrated Sachs- Wolfe effect, a 
fundamental bias between disformally coupled matter 
and baryons and large scale oscillatory features beyond 
the BAO scale. The numerical results and the discussion 
are restricted to the DCDM model, and focus on the role 
of the potential slope 7, which mostly determines the late 
time value of Goff- It remains to be studied whether or 
not similar effects occur in viable models such as the ones 
described in Section [VlB[ and to what extent they might 
be observable by current or future surveys. 

a. Early Dark Energy: Both the baryons and the 
coupled Dark Matter are indistinguishable as long as 
the coupling is negligible. They are equally affected by 
the presence of early Dark Energy (see Figure |4| , which 
produces a departure from the matter domination result 
(5 (X a: EDE increases the expansion rate without cluster- 
ing, reducing the formation of structure. This effect was 
also found for the uncoupled scalar field [351 HO] , and is 
most noticeable for models with higher fiedc (e.g. 7 = 4). 

h. Late Enhanced Growth: The growth of structure 
is enhanced after the transition takes place, consistently 
with the small scale approximation ([72| . The large value 
of Geff/G overcomes the additional friction term, and 
structures form much faster than in the standard CDM 
scenario. Models with less EDE suffer a higher en- 
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Figure 4: Growth of disformally coupled matter and (uncoupled) baryons. Upper panel: Thick purple lines coorespond 
to the disformally coupled matter and thin brown lines to baryons at the same scales. The pressence of early Dark Energy 
affects the slope at early times. Afterwards, the Dark Matter perturbation starts growing rapidly, as deduced in the small 
scale approximation (see Figure [3|. Lower panel: Evolution of the power spectrum at 2; = (thick, solid), 0.3, 0.6, 1 (thick, 
dotted), 2, 4 and 10 (thin, dotted). The gray line corresponds to a reference ACDM model. Units of k are Mpc~^. 



hancement, because the effective gravitational constant 
GeS oc 7^ is larger and the transition occurs earlier (i.e. 
the field takes longer to dominate the energy content). 
The effect from the bump in the effective Newton's con- 
stant associated to the transition is not obvious in the 
evolution of S, and is subdominant with respect to po- 
tential domination. 

The enhanced growth effect is partly canceled by the 
early Dark Energy damping. This degeneracy causes the 
relative resemblance between DCDM power spectra with 
7 = 4 and the DM power for the standard model on 
small scales, but fails anywhere else. It would be worth 
exploring this cancellation in a more systematic way (e.g. 
Markov chain Monte Carlo exploration of the model pa- 
rameters), which would in turn require a better under- 
standing of the baryon bias induced by the coupling (see 
below). However, such an exploration is postponed for 
future work. 

Note that Fourier modes reach the non-linear regime 
earlier due to the enhanced growth. Upon the failure of 
linear perturbation theory, the disformal screening mech- 



anism explored in Section [V] might hide these dramatic 
effects and restore the standard growth, softening the de- 
viations on small scales. Although this seems unlikely to 
save the example model, it might be necessary to take 
the effect into account to obtain a fair comparison with 
observations.^^ 

c. Scale Dependent Growth and Bias: The power 
spectra show scale dependent evolution, as can be seen 
in the different power spectra normalized to the corre- 
sponding ACDM (second line of Figure [5]). In the stan- 
dard model, the linear growth factor is scale independent 
and does not distinguish baryons from Dark Matter. For 
disformally coupled Dark Matter, the scale dependent 
growth follows from the fc-dependent term in the per- 
turbed coupling ( 68 ) . This feature does not appear in 



phenomenological coupled models, in which the growth 



Similar enhanced growth effects have been considered in the con- 
text of quintessence conformally coupled to neutrinos, where the 
necessity of non-linear analysis has been pointed out |99| . 
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Figure 5: Power spectra for disformally coupled matter and (uncoupled) baryons a,t z — for /3 = 2O7 and different values of 
7 (units of k are Mpc~^). First line: Power spectrum for (uncoupled) baryons and (coupled) matter. Second Line: Ratio 
of the power spectra in disformally coupled models relative to ACDM. See the text for a discussion of the different effects. 



of the coupled matter structures is enhanced, but in a 
scale independent way (cf. 

The coupling also modifies the relation between bary- 
onic and Dark Matter structures, since DM couples di- 
rectly to the field while baryons are only indirectly af- 
fected. As baryons are dragged into the potential wells 
created by the coupled matter, they follow a scale depen- 
dent growth pattern, delayed with respect to the domi- 
nant matter component. The resulting bias between the 
two species is shown in Figure |6] The scale dependence 
of the bias vanishes both on super-horizon scales {k/h < 
0.001 Mpc-1) and the small scales {k/h > 0.1 Mpc-1), 
which are well described by the scale-independent Gcs 
(72). The intermediate region shows the interplay be- 



tween the scale dependent growth for the coupled matter 
and the baryons following these structures. 

Since galaxies form out of baryons, this fundamental 
bias modifies the usual DM-galaxy power relation |100j . 
Such a correction needs to be taken into account when 
comparing the observed power spectrum with disformally 
coupled models. On the other hand, other measurements 
of the matter distribution such as weak lensing would 
probe the structures formed by both components, and 



may be used to break the degeneracy. 

d. CMB: Integrated Sachs- Wolfe Effect: The en- 
hancement of the perturbations after the transition 
causes the very large Integrated Sachs- Wolfe effect ap- 
preciated in Figure [6] which becomes most noticeable for 
the models with higher values of 7. The model with 7 = 4 
gives a better fit on the / < 10 multipoles, but departs 
considerably in the range 10 < Z < 200 due to the effect of 
early Dark Energy after recombination. The model with 
less early Dark Energy has the opposite problem: it pro- 
duces a better fit in the intermediate range 100 < / < 200 
due to the lower amount of early Dark Energy, but the 
ISW enhancement explodes at lower multipoles due to 
the higher value of GcS- The different amounts of early 
Dark Energy have an additional effect on the CMB nor- 
malization due to the primary Sachs- Wolfe effect: by re- 
ducing the potential wells that redshift the photons, ficdo 
acts increasing the height of the peaks. 

e. Oscillatory Features beyond the BAO Scale: Os- 
cillatory features can be appreciated in the coupled mat- 
ter power spectrum on very large scales. These are likely 
created as field oscillations on scales near k ^ H{z), 
which are then transferred to the coupled component. 
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Figure 6: CMB power spectrum (left) and bias between baryons and Disformally Coupled DM (DCM) induced by the coupling 
(right). The enhanced growth of Dark Matter structures on small scales produces a very large ISW effect. Note that the 
departures are worse for models with less early Dark Energy (higher 7), as derived in the small scale approximation ( |72[ ). Units 
of k are Mpc~^. 



when the coupling is active. They are most noticeable for 
the models with a large early Dark Energy component, 
e.g. larger field energy density. Although it constitutes a 
distinctive feature of the model, the oscillations are not 
significantly imprinted on the baryonic power spectrum. 
This, together with the large survey volumes necessary 
to explore such scales would make it difficult to detect 
them through LSS surveys. However, the large scale os- 
cillations would be a characteristic signature in models 
where the disformal coupling is universal, in which the 
same effects occur to DM and Baryons. 



B. Viable Scenarios 

The study of cosmological perturbations within the 
Disformally Coupled Dark Matter model (73 75) shows 



very drastic departures in the formation of large scale 
structure, which seem very difficult to reconcile with ob- 
servations. It would be interesting to obtain a more 
precise quantification of these discrepancies through an 
MCMC analysis and explore possible degeneracies (e.g. 
the growth suppression from early Dark Energy and the 
enhancement from the high Gcff). However, it is neces- 
sary to address the existence of alternative, viable sce- 
narios already at this stage. 

Luckily, the action ( 33 ) is very general and there is con- 



siderable room for improvement through different choices 
of the functions A, B and Cc/,. There are at least two pos- 
sibilities 

• Introduce a modulation in the disformal factor 
B{(j)) — >■ f{(j))B{(l)), to make Qq small enough after 
the field enters the slow roll phase. This modifica- 
tion can render 6GeB arbitrarily small, except for 
a relatively short time around the transition (see 
Figure [T]) . This type of models would allow us to 



study the effects imprinted by the transition bump 
without the problems caused by the high Ges at 
late times. 

• Constructing the field Lagrangian using a disfor- 
mal metric, as in the uncoupled model presented 
in Refs. UHl ISO]- In this model the transition to 
slow roll would be partially driven by the scalar 
field Lagrangian itself, and the effects on matter 
may be significantly reduced. In the minimal pre- 
scription, the matter and field Lagrangian are con- 
structed using the same metric ^ and no extra 
parameters are introduced. If this model turned 
out not to be viable, a different disformal metric 
for the field and the coupled matter would offer 
an alternative that is able to interpolate between 
disformal quintessence and the disformally coupled 
Dark Matter presented here (e.g. different disfor- 
mal factors with B^™) = eB^-'t'^). 

Other alternatives could be based on the interplay be- 
tween the conformal and the disformal part of the cou- 
pling. Viable scenarios might be exploited to alleviate 
the claimed problems of ACDM with small scale struc- 
ture formation, such as the tension between Dark Mat- 
ter simulations and observations with regard to both the 
density profiles of Dark Matter halos and for the number 
of predicted substructures inside a given host halo, the 
baryonic Tully-Fisher relation, the constant galactic sur- 
face density or the large scale bulk flows (See reference 
[lOlj for a summary and references therein for further 
details). 

As a final remark, let us note that the enhanced growth 
rushes the modes into entering the non-linear regime at 
earlier times, breaking down the perturbative approach 
followed here. As it was explained in Section |V| the dis- 
formal coupling comes equipped with a screening mech- 
anism, that hides the effects of the additional force on 
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dense environments. Addressing the consequences of the 
disformal screening in a cosmological context would re- 
quire considering the non-hnear backreaction of the field, 
which is not properly captured in the approximations 
considered so far. Chameleon-type theories also show a 
strengthening of the screening when non-linearities are 
properly taken into account [102] . 



VII. DISCUSSION 

The disformal relation provides a generalization of the 
conformal transformation. It has been used to construct 
theories of modified gravity, notably those which produce 
non-trivial effects on null geodesies, such as varying speed 
of light and gravitational alternatives to Dark Matter. It 
also appears in the description of branes embedded in a 
higher dimensional bulk space, in which the scalar fields 
represent the brane position in a certain set of coordi- 
nates. The results of the present work concern the set 
of theories which can be expressed as General Relativity 
plus a matter Lagrangian, which is constructed using the 
disformal metric. This provides a generalization of the 
old school scalar-tensor theories in the Einstein Frame: 
Test particles follow geodesies which explicitly involve 
derivatives of the scalar field, and the energy momentum 
of the field and coupled matter (computed with respect 
to the gravitational metric) is not conserved separately. 

The existence of additional frames, in which only the 
conformal or disformal part enter the matter action ex- 
plicitly, provides novel connections between scalar-tensor 
theories of gravity. In particular, it is possible to restore 
the theory to a Jordan Frame representation by revers- 



ing the disformal relation, as was shown in Section [ill A 



By doing so, the transformed Einstein-Hilbert term is 
shown to be equivalent to a quartic DBI Galileon La- 
grangian when expressed in a frame in which the disfor- 
mal coupling is pushed towards the gravitational sector. 
The resulting theory has the correct Horndenski form 
([T]), ensuring the second order nature of the equations 
and the avoidance of Ostrogradski's ghosts. In particu- 
lar, it introduces a derivative coupling between the scalar 
field and gravity, together with higher derivative self in- 
teractions. These endow the theory with the Vainshtein 
screening mechanism, which allows the field to cause ef- 
fects on cosmological scales while remaining undetectable 
in the Solar System. 

The equivalence between certain higher derivative the- 
ories (such as DBI Galileons with conformal or mini- 
mal coupling to matter) and disformally coupled the- 
ories with an Einstein-Hilbert gravitational sector pro- 
vides new means to analyze this type of models. Al- 
though the equations for disformally coupled theories 
are rather involved, they are much simpler than higher 
derivative Horndenski theories in the Jordan Frame and 
highlight the properties of the different terms. Hence 
the analysis of disformal theories in the Einstein Frame 
can be regarded as equivalent to (at least) some scalar- 



tensor theories featuring the Vainshtein screening mecha- 
nism. The kinetic mixing between the coupled degrees of 
freedom makes it necessary to solve for the second time 
derivatives of the different components. Although this 
generally requires some assumptions, e.g. about the met- 
ric, a general equation without derivatives of the energy- 
momentum tensor can be obtained (45 1. Once solutions 



are computed, it is possible to restore to the Jordan met- 
ric in order to interpret the results and compare to ob- 
servations. 

In high density environments (as measured by the con- 
dition Bp ^ 1) the field does not feel the presence of dis- 
formally coupled, non-relativistic matter. This provides 
a novel disformal screening mechanism, which is distinct 
from screening mechanisms based on the field potential 
(Ghameleon and Symmetron jl03) ). Our mechanism re- 
lies on the existence of a well defined limit p — > c» in the 
scalar field equation, -'^'^ given by equation (54), for which 



the field evolution is independent of the matter distribu- 
tion and the field gradients (up to effects of order ~ p/ p, 
v/c). If the conformal part A is negligible, only a fric- 
tion term remains and the field coupling density ( 43 1 is 
a decreasing function of time. As it evolves below its 
cosmological value (provided V' <Q and B' /B > 0), the 
effects of the coupling are suppressed by a factor ~ p^ / p 
and the theory is consistent with precision gravity tests. 
Potentially detectable signatures may be obtained in the 
presence of matter velocity flows, radiation pressure or 
relativistic matter, strong gravitational fields or gradi- 
ents of cosmological origin. 

The disformal screening mechanism may be related 
to the Vainhstein effect, which suppresses the gradients 
of the scalar field and hides the additional force near 
massive sources due to the higher order derivative self- 
interactions. In the disformal case, the screening relies 
on the kinetic mixing between the scalar field and the 
coupled degrees of freedom. This ultimately allows the 
existence of a well defined p — ^ oo, non-relativistic limit 
in which the field is free, or only subject to friction (up 
to conformal interactions). Therefore, both the disformal 
and Vainshtein mechanisms belong to the kinetic screen- 
ing category, as they rely on the form in which the field 
derivatives occur in the action. Moreover, the two cases 
might indeed be the same physical effect expressed in dif- 
ferent frames, due to the equivalence between the disfor- 
mally coupled theories and a particular higher derivative 
sector of the Horndenski Lagrangian. 

The disformal coupling offers interesting possibilities 
to build models for cosmic acceleration. In the FRW 
approximation, the same properties that gave rise to 
the disformal screening mechanism make the background 
coupling approximately proportional to the Dark Energy 



One may show that this limit is independent of the assumption 
of canonical kinetic term for the scalar field wo adopted for most 
of our discussion in this article. 
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density (63) rather than to the coupled matter energy Frame pressure if Bp > 1, as it was discussed in Section 



density. This provides a concrete reahzation of a class of 
interacting Dark Matter models which have been exten- 
sively studied using phenomenological parameterizations. 
The equations for linear perturbations around FRW con- 
tain scale dependent terms. These are absent in the pure 
conformally coupled case, and have hence the potential 
to distinguish the two possibilities. An analytic equation 
for the coupled matter perturbations was derived in the 
small scale limit. On top of an additional friction term, 
the effect of the fifth force can be encapsulated in the def- 
inition of an effective gravitational constant ( 72 1 which 
depends on the background coupling factor. 

In order to investigate the cosmological implications 
of a disformal coupling in a simple setting, a Disformally 
Coupled Dark Matter (DCDM) example model was pro- 
posed. This has the advantage of avoiding the subtleties 
of the Einstein Frame description, since gravity, baryons 
and photons share the same physical metric. A DCDM 
model with exponential functions and no conformal cou- 
pling ( 73]|75 1 provides a Dark Energy model that tracks 
the dominant energy component at early times. When 
the coupling to Dark Matter becomes active, the scalar 
field enters a slow roll phase in order to dynamically avoid 
a singularity of the disformal metric. The free parame- 
ters can be constrained by observations, and the model is 
successful at the background level. When perturbations 
are included, the DCDM model introduces a series of 
new effects. The effective gravitational constant for this 
model is too large, due to the persistence of the coupling 
at late times and the domination of the scalar field in the 
energy budget. This causes a too large enhancement of 
the growth factor, which affects the normalization of the 
DM and baryon power spectra, producing a very large 
ISW effect. Scale dependent effects are reflected on mat- 
ter oscillatory features on very large scales and a scale 
dependent bias between the coupled Dark Matter and 
(uncoupled) baryonic component. 

There is considerable freedom in the model to produce 
cosmologically viable scenarios. Models of the DCDM 
type with less dramatic growth of perturbations can be 
constructed by modifying the functional dependence of 
the disformal coupling (e.g. tuning it to become neg- 
ligible after the transition to slow roll), the scalar field 
Lagrangian (e.g. constructing it with a disformal metric), 
or perhaps by the interplay between the conformal and 
the disformal parts of the coupling. Unlike in the con- 
formal case, the disformal coupling affects ultrarelativis- 
tic species and hence variations or extensions of DCDM 
may postulate or include disformally coupled neutrinos 
or photons. 

Another phenomenological direction is to consider the 
disformal screening mechanism in detail. The results 
presented here considered a purely disformal coupling, 
monotonically increasing with the field. Certainly, in- 
cluding a conformal coupling and more general functional 
forms is of interest. These considerations might help 
to avoid the gradient instability caused by the Einstein 



IV A Once different set-ups are formulated, it is worth to 
explore the observable signatures for the model by quan- 
tifying the effects outhned at the end of Section [VA| 

The dependence of the free functions in the Horndenski 
Lagrangian Q [5| on the field kinetic term X has a very 
special role, as it relates the coupling to gravity to the 
coefficients of the second derivative field terms. There- 
fore, it would be worth considering the transformations 
between frames in the more general case in which the 
disformal relations are allowed to depend on X. The 
computation of the Ricci scalar associated to this gen- 
eral disformal metric would provide the Jordan Frame 
representation of the most general scalar-tensor theory 
that accepts an Einstein Frame description. Since the 
equations simplify considerably in this frame, the phe- 
nomenology of these theories would be relatively easy to 
address. 

Finally, the existence of a well behaved p — >■ oo classical 
limit in the field equation suggests that disformally cou- 
pled theories might introduce new interesting features for 
the physics of gravitational singularities and other high 
energy regimes. The implications of kinetic mixing for 
the formation of black holes or the origin of the universe 
is beyond the scope of the present work, but it might pro- 
vide a fruitful exploration to pursue in the future. This 
discussion provided just a glimpse to the potential appli- 
cations of the disformal relation. As a generalization of 
the conformal case, which was very central to the devel- 
opment of gravitation and cosmology in the 20*^ Century, 
the use of disformal transformation might provide novel 
ways to address the gravitational physics of the 21*'^ Cen- 
tury. 
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Appendix A: Disformal Relations 

Consider the disformal relation between two metrics, 
specified by the two scalar functions A, B, and a vector 



g^u = Agf_,^ + Bb^,K ■ 



(Al) 



20 



The inverse metric can be found by contraction 
1 



9 



A 



5 



where 



B 



7 



(A2) 



(A3) 



where 



D 



_ 5r'' _ 1 



7^" A 



(A8) 



The inverse relations are provided below for complete- 
ness 



and b^ = g^'^b^Jji, = W'b^. The determinant of the barred 
and unbarred metrics are related 



A, 



lAB 



A\ I 



B 



62, 



(A4) 



The above relation is derived in Appendix C of Ref. [3T] . 

It is possible to write the relation of stress energy mo- 
mentum tensor (associated to a Lagrangian ^J—gC) in 
the two metrics by using the chain rule 



\J-g Sgf,^ V 9 S9tJ.iy 



2 Sj^C) 

^ 5gai3 



(A5) 

By identifying the quantity in brackets as T^^'^ and using 
( A4 ) , the following relation follows 



B 



A-Bb^ 



(A9) 



(AlO) 



(All) 

where h^' = g'^'^K. Note that = b^, bt" = B/{Aj^)bi' 
and 7^62 = Bb^. 



B 

A^ 
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(A6) 



The equivalent relation with lower indices is considerably 
more involved 



(A7) 



Appendix B: Disformal Geodesies 



The expression for the disformal connection (11) can 



be expanded in terms of the functions in the disformal 
metric 



Kp + 5l^ogA^,)--\ogA'^g^, 



B 



A 



2 r I f I 

A(a<^,/3) - 7,<l^'^A,xgo.p ~ 2X B - -(j}'^B^x(l),o.(l),p 



(Bl) 



A-2BX s-rises from the inverse barred metric, eq. (A2|). The first term is just the connection of the 
unbarred metric, and the two following terms arise from the purer 



Here ^2 — §_ 

ly conformal transformation involving derivatives of 
A. The fourth term and the second line contain the first order derivative terms from the disformal contribution to 
the metric B. The third line shows the second order derivative terms VV0. 



Appendix C: General Perturbations 



The coupling density perturbation that enters the linear equations (67 65) in the case where both the conformal 



and the disformal parts of the coupling are relevant has the rather complicated form 



6Q = QpS^c + Q4>H + Qd<t,6(t) + Q** + , 



(CI) 
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where 



B •,\ A' (A - 2502) + __ M2BV' + QBHcj)) 

'.P - ^ll-T'^") 7 — 

^ 2(A + Bp- B<l)^ 



Q4> = ( (1 - 2^(/.2)A" + _ 25 



V' 



A'' 

'A 



Qd<j, 



2A{A + Bp- B(/)2) 

(A^ - 2B{2A + Bp)(j)^ + 25^04) + A' (^-2B{V' + iH4>) + B' {p + 20^ 

2AB'{V' + 3ff^) + A2S'2^2(^ _ ^2)\ ^ , 

/ 2{A + Bp - B<j?Y 

B{A + 2Bp)^(j> - iA + Bp)B'^ + B (2BV'(j> + 3H{A + Bp 

A + Bp- Be 



B{A + 2Bp)^(j) -{A + Bp)B'4) + 2B {?,H{A + Bp) + BV'c, 
{A + Bp- B4>^)^ 

3Bp(}) 



(C2) 



(C3) 



(C4) 



(C5) 



(C6) 



A + Bp-B<p^ 

The equations for the perturbations can also be found in Ref. [90| for the case in which the coupled fluid is allowed 
to have pressure. 



Appendix D: Dynamical system analysis 

It is useful to reformulate the system in terms of the dimensionless variables 



8nGp _ UttG d(t) _ %ttGV _ BH^ 

and using the e-folding time N = log a as the time variable. The Friedmann constraint then reads 

l = n + X^ + Y . 
For concreteness, we assume the exponential forms 

A = Aoe""^/*^" , B ^ Boe^"^/^^" , V = Voe''^'''/'^''' , 
with Mp = (47rG'/3) . Note that when a ~ p one has a simpler form 



(Dl) 

(D2) 
(D3) 



1 

M4' 



We use Eq. ( D2 1 to eliminate f2 from the system. Equations of motion can then be rewritten in terms of the 

(D4) 



remaining variables as 

A' = aAX , 

X' 



1 r 



A-3Z {3X^ + Y -l) A {6X^ + X^a - &X{Y + I) + {Y - l)a - 2Yj) 
- 6XZ {9X* + X^{2a - P) - QX^{Y + 1) + X{Y - l){2a - (3) - 2XY-i - 'i{Y - 1) 

Y' = 3 ( 1 + + ^7^- y ) y , 



Z' = -^[l + X^ - -(3X-Y]Z . 



(D5) 
(D6) 



The fixed points are: 
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• Matter domination: O = 1, X = 0. This solution is always a saddle point, since the eigenvalues corresponding 

to it arc (3, —3, —ajl). 

• Scaling solution: f2 = 1 — X = — a/2. This point in the phase space is never an attractor. 

• Conformal scaling solution: f2 = 2 { ^^^^^2^^'^ ) > ^ = ■ia-^i ' "^^^ eigenvalues are 

6(^ + 7) ± v^(a - 27)2 (- (4a2 + 63) 7^ + a {a^ - 54) 7 + 45^2 + Aar' + 1296) - 3(a - 27)(a - 7) 
^ a -27 ' 2(q! - 27)2 ^" 

The general form of the stability condition is too messy to write down here. 

• Kinetic domination: f2 = 0, X = ±1. This solution is stable given ±a < —6 A ±/? < 6 A ±7 < —6. 

• Scalar dominated solution: = 0, X = —7/6. The stability conditions for this fixed point are modified in the 
presence of the disformal coupling. To be explicit, this solution exists and is stable if either a+y > 27 A7 < 6A 

> OA7 > 0)V(-6 < /3 < OA/3+7 > 0)) or 0;+^ < 27A7+6 > 0A((7 < OA/3 < 0)V(/3 > OA^+7 < OA^ < 6)). 

• Disformal scaling solution: O = (36 - /32 ± ^^^"^ - m)l\%, X = {/3 ± ^//S'^ - 36)/6. This fixed point exists 
when > 6, but it is never stable. 

I 

The last one of these is new and it exists when A = 1. 
The plus-branch is physical, given /3 > 6. 
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